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ABSTRACT 

In  situ  measurements  of  ion  charge  states  can  provide  unique  insight  into  the 
heating  and  evolution  of  coronal  mass  ejections  when  tested  against  realistic  non¬ 
equilibrium  ionization  modeling.  In  this  work  we  investigate  the  representation 
of  the  CME  magnetic  field  as  an  expanding  spheromak  configuration,  where 
the  plasma  heating  is  prescribed  by  the  choice  of  anomalous  resistivity  and  the 
spheromak  dynamics.  We  chose  as  a  test  case,  the  19  May  2007  CME  observed 
by  STEREO  and  ACE.  The  spheromak  is  an  appealing  physical  model,  because 
the  location  and  degree  of  heating  is  fixed  by  the  choice  of  anomalous  resistivity 
and  the  spheromak  expansion  rate  which  we  constrain  with  observations.  This 
model  can  provide  the  heating  required  between  l.li?Q  and  earth  orbit  to  produce 
charge  states  observed  in  the  CME  flux  rope. 

However  this  source  of  heating  in  the  spheromak  alone  has  difficulty  account¬ 
ing  for  the  rapid  heating  to  Fe8-11+  at  lower  heights,  as  observed  in  STEREO 
EUVI  due  to  the  rapid  radiative  cooling  that  occurs  at  the  high  densities  involved. 
Episodes  of  heating  and  cooling  clearly  unrelated  to  spheromak  expansion  are  ob¬ 
served  prior  to  the  eruption,  and  presumably  still  play  a  role  during  the  eruption 
itself.  Spheromak  heating  is  also  not  capable  of  reproducing  the  high  Fe  charge 
states  (Fe16+  and  higher)  seen  in  situ  exterior  to  the  flux  rope  in  this  CME.  Thus 
while  the  spheromak  configuration  may  be  a  valid  model  for  the  magnetic  topol¬ 
ogy,  other  means  of  energization  are  still  required  to  provide  much  of  the  rapid 
heating  observed. 


1.  Introduction 

The  eruption  mechanism  of  coronal  mass  ejections  (CMEs)  is  currently  an  extremely 
active  area  of  research.  A  variety  of  mechanisms  have  been  discussed  in  the  literature. 
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Most  of  these  require  some  form  of  magnetic  reconnection  during  and  after  the  eruption. 
The  associated  thermal  heating  should  be  an  important  diagnostic  of  the  eruptive  process, 
including  the  eruption  dynamics  and  ejecta  geometry,  but  until  now  has  been  relatively 
unexploited.  Such  heating  gives  rise  to  UV-X-ray  emissions  during  the  eruption,  that  can 
be  remotely  sensed,  and  also  leaves  its  imprint  in  the  charge  states  of  ions  detected  in  situ 
near  1  AU,  as  their  ionization  state  responds  to  the  newly  heated  plasma.  These  techniques 
give  complementary  insights  into  the  heating  and  evolution.  Spectroscopy  is  in  some  sense 
more  direct,  being  a  prompt  signature  of  the  eruption,  but  with  existing  instrumentation 
has  only  been  possible  relatively  close  to  the  solar  surface  where  the  signal  is  strong.  Ion 
charge  states,  on  the  other  hand,  are  routinely  collected  and  CMEs,  and  their  constituent 
parts  (forward  shock,  flux  rope,  etc)  can  be  easily  identified.  Further  interpretation  of  such 
data  require  modeling  of  the  non-equilibrium  ionization  (NEI)  balance  as  the  CME  plasma 
expands  outward  into  the  ambient  solar  wind.  In  this  paper  we  explore  the  consequences  of 
a  different  source  of  thermal  heating,  that  due  to  anomalous  resistivity  within  a  spheromak 
solution  for  the  CME  geometry  as  it  expands  into  the  solar  wind.  We  compare  the  predictions 
of  such  a  model  to  STEREO  and  ACE  observations  of  the  19  May  2007  CME.  The  structure 
of  the  paper  is  as  follows.  We  describe  the  status  of  our  understanding  of  ion  charge  states 
observed  in  CMEs  in  more  detail  in  section  2.  Section  3  describes  the  2007  May  19  CME, 
and  some  of  our  motivation  for  exploring  a  spheromak  model.  Section  4  gives  the  theory 
of  spheromak  resistive  expansion,  with  some  discussion  of  just  how  closely  the  spheromak 
model  used  here  may  represent  a  real  CME.  Section  5  outlines  the  NEI  model  incorporating 
the  spheromak  resistive  heating,  and  we  conclude  in  section  6  with  a  summary  of  our  results 
in  the  context  of  previous  research  and  discuss  future  modeling  improvements. 


2.  NEI  modeling  of  CME  charge  states 

Interplanetary  CMEs  or  ICMEs  have  many  well-known  signatures.  Of  these,  in  situ 
measurements  of  ion  charge  state  distributions  hold  a  unique  potential  for  diagnosing  the 
conditions  throughout  the  eruption.  The  charge  states  evolve  through  ionizations  and  re¬ 
combinations  as  a  function  of  temperature  and  density  up  to  heights  of  3  to  6  solar  radii 
heliocentric  distance  but  are  “frozen  in”  thereafter.  Since  low-Z  ions  “freeze-in”  at  higher 
densities  (i.e.  smaller  heights  or  radii)  than  high-Z  ions,  they  provide  complimentary  con¬ 
straints  on  the  temperature  evolution. 

A  number  of  studies  have  examined  the  existing  solar  wind  composition  data  and  made 
inferences  on  freeze-in  temperatures  based  on  computed  ionization  distributions  appropriate 
to  coronal  equilibrium  (e.g.  Zurbuchen  and  Richardson  2006,  Zurbuchen  2004,  Lepri  and 
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Zurbuchen  2004,  Lepri  et  al.  2001,  Henke  et  al.  2001,  Gloeckler  et  al.  1998,  Henke  et  al. 
1998).  However  Rakowski,  Laming  &  Lepri  (2007)  were  the  first  to  conduct  time-dependent 
modeling  of  the  ion  charge  state  distributions  of  various  elements  detected  in  situ  in  order 
to  draw  quantitative  conclusions  regarding  the  thermal  input  and  initial  conditions  in  the 
corona  during  the  CME  eruption.  Rakowski.  Laming  fe  Lepri  (120071  modeled  the  charge 
states  of  8  ICMEs  from  the  event  list  of  Lynch  et  al.  (120031),  chosen  to  cover  a  range  of 
velocities. 


Rakowski,  Laming  fe  Lepri  (20071)  considered  a  simplified  geometry  for  the  CME  ejecta 


and  parameterized  the  evolution.  The  model  velocity  and  expansion  were  based  on  the  ob- 


evolution  consists  of  three  phases:  initiation,  acceleration  and  propagation 

(Zha 

ng  et  al. 

2001, 

2004; 

Zhang  &  Dere 

2006; 

Sui  &  Holman 

2003; 

Lin  et  al. 

2005; 

Sui  et  al. 

2005 

).  Heat- 

ing  of  the  CME  core  plasma  during  the  acceleration  phase  was  explored  assuming  a  heating 
rate  for  the  CME  plasma  proportional  to  the  rate  of  kinetic  energy  increase,  i.e.  a  constant 
fraction  QE/KE  during  the  acceleration  up  to  a  final  velocity  Vf.  In  5  out  of  8  sample  CMEs 
studied,  the  dominant  Fe  charge  states  were  neon-like  (16+)  or  higher,  indicating  that  high 
temperatures,  comparable  to  flares  (  107K),  are  involved.  Starting  the  plasma  from  this 

temperature  and  allowing  the  ions  to  recombine  as  they  expand  could  often  account  for  the 
Fe  ionization  balance,  with  peaks  around  Fe16+  and  Fe8+  (the  Ne-like  and  Ar-like  charge 
states,  which  have  small  recombination  rates  to  the  next  charge  states  down,  hence  popu¬ 
lation  “bottlenecks”  here).  However  the  lower-Z  elements  placed  a  limit  on  the  maximum 
starting  temperature  (at  least  if  assuming  ionization  equilibrium  in  the  seed  plasma).  Above 
~  2.5  x  106K,  O  would  have  been  mainly  08+  instead  of  07+  and  06+  as  observed  and 
would  not  recombine  significantly  during  the  CME  evolution.  Evidently  plasma  must  start 
out  much  cooler,  and  be  further  heated  as  the  CME  accelerates.  Thermal  energy  inputs  on 
the  order  of  2  to  10  times  the  kinetic  energy  were  needed  in  7  out  of  8  events  studied. 


Similar  conclusions  to iRakowski.  Laming  fe  Lepri  (120071)  about  the  thermal  energy  in- 


put  to  CMEs 


lave  been  reached  from  analysis  of  ultra-violet  spectra  taken  by  SOHO/UVCS. 

a  thermal 


Akmal  et  al.l  (200ll)  studied  a  480  km  s  1  CME  observed  on  1999  April  23,  and  fine 


energy  comparable  to  the  bulk  kinetic  energy  of  the  plasma.  Ciarayella  et  ah  (20011)  give 
similar  results  for  the  260  km  s~x  1997  December  12  CME.  More  dramatically,  Lee  et  al. 


(120091)  studied  the  2001  December  13  event,  and  using  the  combination  of  [O  V]  density  di 
agnostics  and  non-equilibrium  ionization  modeling  of  O  VI  found  that  75%  of  the  magnetic 


energy  must  go  into  heat  to  match  the  UVCS  observations.  IRakowski.  Laming  fc  Lepri 


(2007)  found  a  heat  to  kinetic  energy  ratio  of  at  least  6  for  the  same  event. 


The  motivations  behind  the  current  work  are  twofold.  First,  to  test  if  our  non-equilibrium 
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ionization  modeling  of  the  charge  state  distributions  can  hold  up  to  the  scrutiny  of  additional 
constraints  on  the  eruption  dynamics  coming  from  STEREO  observations  in  the  2007  May 
19  CME,  as  well  as  observations  regarding  the  charge  states  closer  to  the  sun.  Secondly  we 
explore  whether  a  more  physically  and  theoretically  specified  heating  model  when  incorpo¬ 
rated  into  our  non-equilibrium  modeling  can  reproduce  the  charge  states  of  a  real  event  with 
reasonable  physical  parameters. 


3.  The  2007  May  19  CME 


The  2007  May  19  CME  studied  in  this  paper  is  qualitatively  different  from  those  CMEs 


discussed  above.  As  one  of  the  first  major  CME  eruptions  in  the  STEREO  era  t 
May  2007  event  has  been  intensely  studied  by  multiple  authors 


re  19 


Veronig  et  ah 

2008; 

Li  et  ah 

2008 

Liu  et  ah 

2008; 

Bone  et  ah 

2009; 

Kerdraon  et  ah 

201oj).  It  began 

(Gooa 

swamv  et  ah 

2009 

2009; 

Liewer  et  ah 

2009 

It  began  as  a  filament  eruption  from  the  active 
region  AR10956  and  was  detected  as  an  ICME  by  STEREO  B,  ACE  and  possibly  STEREO 


filament  which  eventually  erupted  represents  6%  of  the  total  unsigned  flux  change  during 
this  two  day  period. 

Our  model  of  the  heating  during  the  CME  eruption  and  evolution  through  interplanetary 
space  will  need  to  match  both  the  charge  states  seen  in  the  filament  during  eruption  and 


detected  in  situ,  as  well  as  the  density  and  speed  of  the  ejected  material.  Liewer  et  ah 


(120091)  present  the  STEREO  EUV  images  showing  the  heating  of  the  filament  in  the  hours 
and  minutes  immediately  preceding  the  eruption.  Hot  spots  are  seen  in  the  EUVI  17lA 
filter  (mainly  Fe  IX  and  XII)  at  heights  near  1.07  R0.  The  magnetic  cloud  is  determined 
to  have  mainly  interacted  with  STEREO  B,  which  penetrated  the  center  during  most  of 
May  22  (DoY  142),  while  STEREO  A  only  passed  through  the  magnetic  cloud  periphery 
(Kilpua  et  ah  20091.  The  in  situ  Fe  charge  states  seen  with  PLASTIC  on  STEREO  B  and 
A  are  shown  in  Figure  1 0.  Charge  state  data  from  ACE  for  multiple  ions  is  shown  in  Figure 
2.  The  highest  Fe  charge  states,  Fe14+,  Fe15+,  are  only  seen  in  ACE  and  STEREO  A  which 
were  not  centrally  located  in  the  event,  but  rather  on  its  flanks.  Furthermore  the  appearance 


1Links  to  the  level  2  data  are  available  at  http://fiji.sr.unh.edu/ 
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of  the  highest  charge  states  is  after  the  cloud  passage  as  defined  in  Liu  et  ah  020081). 


The  shock  front  around  the  2007  May  19  CME  launched  at  958  km  s~x,  the  average 
speed  in  transit  to  1  AU  was  around  700  km  sy1  and  the _speed  of  the  MC  as  it  passed 
the  STEREO  B  spacecraft  was  482  km  s-1  fKilpua  et  al.  2009).  The  filament  that  erupted 
from  the  active  region  at  about  the  same  time  only  rose  at  an  average  speed  of  103  km  s^1 
(jLiewer  et  al.l  2009),  but  it  could  have  been  dragged  out  by  the  CME  and  the  fast  solar 
wind  to  reach  a  similar  coasting  velocity  as  the  ICME.  It  had  no  distinct  initiation  or 
acceleration  phase,  and  clearly  does  not  fit  within  the  phenomenology  described  above  and 


m 


Rakowski,  Laming  &  Lepri  (120071) .  which  leads  us  to  consider  a  spheromak  model  for  the 


CME  evolution. 


4.  Resistive  Spheromak  model  of  CMEs 
4.1.  Spheromak  Field  Configurations 


Spheromaks  are  well  known  force-free  configurations  of  plasma  satisfying  condition 
V  x  B  =  oB  with  spatially  constant  a.  They  are  solutions  of  the  Grad-Shafranov  equa- 
tion  in  spherical  coordinates  with  the  poloidal  current  being  a  linear  function  of  the  flux 
Chandrasekhar  fe  Kendall  1957).  The  basic  spheromak  solution  is 


Br  = 
Be  = 

Btf,  = 


~B< 


ar 

j i  (ar)  +  arj[  (ar) 


ar 

Boji  (ar)  sin  6 


sin  6 


(1) 


where  the  spherical  Bessel  function  j\  can  be  expressed  in  terms  of  elementary  functions, 
ji(x)  =  sinx/x2  —  cosx/x.  The  parameter  a  is  related  to  the  size  of  the  spheromak  R, 
defined  by  the  surface  where  radial  magnetic  field  is  zero,  given  by  the  solution  of  j\  =  0; 
R  =  Ca/a  with  Ca  =  4.49. 


The  spheromak  model  of  CME  topology  is  an  attractive  alternative  to  the  usual  flux 
rope  model.  It  is  self-contained  with  no  unknown  “length  of  flux  rope”  parameter  to  main¬ 
tain  connectivity  to  the  sun  and  the  heating  within  it  is  mathematically  specified  once  the 
anomolous  resistivity  is  chosen.  The  torus  structure,  seen  on  edge,  may  even  resemble  the 
filament  shape  in  some  CMEs  (IKataoka  et  al.l  120091) .  Previously,  iLvutikov  fc  Oourgouliatos 
(2010)  found  ideal  self-similar  solutions  for  expanding  spheromak,  with  electric  fields 


r  a 

E  =  — er  x  B 

c  a 


(2) 
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(dot  denotes  differentiation  with  respect  to  time)  and  the  corresponding  non-radial  velocity 


field 


(3) 


In  this  Section  we  discuss  continuous  heating  of  an  expanding  spheromak,  generalizing  the 
previous  analysis  to  non-ideal  self-similar  expansion  with  E  •  B  ^  0. 

Physically,  the  parallel  component  of  the  electric  field  should  be  related  to  the  cur¬ 
rent  density  through  Ohm’s  law.  Formally,  the  procedure  described  below  breaks  down  the 
assumption  of  self-similarity,  since  the  value  of  the  parallel  electric  field  component  is  not 
linearly  proportional  to  the  current  density.  Still,  we  assume  that  resistivity  plays  a  sub¬ 
dominant  role,  so  that  the  expansion  remains  approximately  self-similar.  In  other  words,  we 
assume  that  resistivity  leads  to  small  deviation  from  the  ideal  self-similar  expansion.  The 
resulting  dynamics  remain  self-similar,  by  assumption,  but  are  a  little  different  from  the 
ideal  case.  We  are  interested  not  in  the  detailed  properties  of  local  resistive  heating,  but 
in  general  scaling  relations.  As  we  show  below,  the  requirement  that  an  expanding  CME 
dissipates  some  of  its  initial  energy  in  the  form  of  heat,  detectable  by  its  effect  on  element 
charge  states  by  the  time  it  reaches  Earth  orbit,  can  be  used  to  estimate  the  anomalous 
resistivity. 

Expressing  electric  field  in  radiation  gauge  E  =  —dtA /c,  where  A  is  vector  potential, 
we  find  from  the  force-free  condition, 


dr(rAe)  -  deAr  =  arA 


(4) 


This  equations  highlights  two  important  points.  First,  it  shows  how  the  time-dependences 
of  the  poloidal  components  of  the  vector  potential  are  related  to  the  toroidal  components, 
Ag,  Ar  oc  a(t)A^  (this  is  related  to  the  conservation  of  poloidal  and  toroidal  fluxes).  Second, 
we  get  one  equation  for  the  two  functions  Ar,  Ag.  Since  spheromak  solutions  are  linear,  a 
general  solution  is  a  linear  combination  of  the  two  solutions.  The  two  particular  solutions, 
which  in  a  static  case  correspond  to  two  equivalent  choices  of  the  poloidal  components  of 
vector  potential,  in  a  time  varying  case  correspond  to  different  physical  processes.  Accord¬ 
ingly,  there  are  two  types  of  solutions,  corresponding  to  two  choices  of  vector  potential  in 
Eq.  (|4]),  Ar  7^  0  and  Ag  ^  0.  We  consider  them  in  turn. 
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4.2.  Resistive  expansion  of  first  type 


Let  us  assume  that  B $  scale  with  time  as  an  arbitrary  function  of  the  expansion  param¬ 
eter,  B<p  oc  f(a/a0).  The  first  type  of  solutions  has  Ar  =  0  and,  from  equation  1, 


Aq  —  —B0f(a/ao) 


(sin  {ar)  -  far  ^ dz ) 


a2r 


sin  6. 


(5) 


The  first  electromagnetic  invariant  E  •  B  oc  (aj  —  2 d/) .  Thus,  if  /  =  (a/oo)2  the  flow  is 

\  "  7 

ideal,  consistent  with  the  scaling  chosen  in  Lyutikov  &  Gourgonliatosl  (2010).  More  generally, 
if  we  allow  /  =  [a/ a o)2+m,  we  find  the  parallel  electric  field 


m2  By(a  /  a0) 


2m 


a 2  sin4  6 
ao°2 


Fi(ra,  6) 


(6) 


where  T\(ra,  6)  is  a  lengthy  function  not  given  explicitly  here.  The  case  of  m  =  0  corresponds 
to  ideal  expansion.  Since  we  expect  that  reconnection  eliminates  magnetic  flux,  physically 
realizable  solutions  correspond  to  m  >  0.  Magnetic  flux  decreases  as  am  oc  R~m. 


4.3.  Resistive  expansion  of  second  type 

The  second  type  of  dissipative  expansion  corresponds  to  Ae  =  0  and 

Ar  =  B0rf(a/a0)j1  cosd 
The  parallel  electric  field  in  this  case 

rv ^ 

El  2  =  Bl(a/a0f^—F2(ra,e). 

CYqC 


(7) 


(8) 


For  this  solution  the  resulting  electric  field  has  a  component  along  magnetic  field  for  any 
m.  The  resistive  effects  in  these  solutions  are  induced  exclusively  by  spheromak  expansion. 
Note,  that  this  is  different  from  resistive  decay  of  a  stationary  spheromak,  which  proceeds 
homologously.  Thus,  expansion  of  the  spheromak  generally  leads  to  appearances  of  parallel 
electric  fields  and  associated  dissipation,  which  can  proceed  much  faster  than  resistive  decay 
of  stationary  spheromaks. 


4.4.  Expansion-driven  dissipation 

Let  us  concentrate  on  the  resistive  expansion  of  first  type,  which  produces  equatorial 
dissipation  concentrated  in  the  flux  rope,  as  in  the  CME  observations  we  attempt  to  model. 


The  dynamical  models  described  above  relate  the  resistivity  to  the  overall  dynamics.  The 
volume  integral  over  E2  (see  Eq.  (EJ))  is 

E\dV  =  0.033 m2B2R3  —  f  J  (9) 

which  leads  to  the  dissipation  rate 

/E^1  f  /-?  \  2?77-— (—  1 

-J-dV  =  0.033 mBlBl  (-^)  |  (10) 

with  the  identification  of  the  resistivity  r]  =  mRv/c 2 .  This  is  justified  by  integrating  the 
energy  dissipation  rate  above  between  R0  and  R  to  give  a  total  dissipated  energy 

-dr  =  0.033 —^vRl  (l  -  ( \  ~  0.033— x  2m  In  (11) 

v  2rj  c2  0  ^  \RJ  )  2 r]  c2  0  \R0J  y  J 

where  the  last  step  assumes  ( R  —  R0)  / R0  «  1.  With  total  spheromak  (magnetic)  energy 
U  =  1.57  x  10~2BqRq/ R  x  (R0/R)2m,  the  difference  in  energy  between  R0  and  R  is  5U  = 
U  (R0)  (l  —  (R0/ R)2m+1) ,  of  which  a  fraction  2m/  (2m  +  1)  is  dissipated.  Writing 

5U  ~  2 mU  (R0)  In  (  ■ E  )  =  0.033 B2R30  In  (  E  )  (12) 

\Ro  J  \Ro  J 

and  equating  with  ^ dr  in  equation  11  yields  rj  =  mR^v/c2.  This  is  derived  assuming  an 
infinitesimal  expansion.  The  generalization  to  arbitrary  expansions  gives 

r)  ~  mRv/c2.  (13) 


In  the  second  case  the  dissipation  rate 

/  p  \  2(m+l) 

E  •  J dV  =  B2R30  VF2(m )  (14) 

where  F2(m)  ~  0.033m  —  5  x  10-4,  very  similar  to  the  first  case  when  integrated  over  the 
spheromak  volume.  The  distribution  of  the  heating  is  different  in  the  two  cases,  as  illustrated 
in  Figure  [3j  The  amount  of  energy  available  for  heating  therefore  depends  on  the  value  of 
the  anomalous  resistivity  and  varies  with  position  in  the  spheromak  about  an  average  value. 
The  distribution  of  volumetric  heating  rates  for  the  first  type  of  dissipative  solutions  (of 
most  interest  here;  see  below)  is  given  in  Figure  IU 

To  what  extent  can  a  spheromak  approximate  a  real  CME?  Flux  rope  structures,  similar 
to  the  equatorial  portion  of  a  spheromak  are  frequently  observed  in  situ,  through  the  rotation 
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of  the  magnetic  field  as  the  spacecraft  penetrates  the  CME  plasma.  The  precise  origin  of  the 


flux  rope  is  not  clear, 
the  photosphere  (e.g.. 


Many  authors  assume  that  such  a  structure  emerges  fully  formed  from 


Fan  fc  Gibson 


breakout,  (e.g.  lAntiochos  et  ah 


20041  iMagara  fc  Longconel  120031).  Alternatively  in  the 


19991:  iLvnch  et  al.ll2004l).  tether-cutting  (e.g.  iMoore  et  al. 


20011  or  flux  cancelation  (e.g.  Linker  et  ah  2003h  models  the  flux  rope  forms  as  the  CME 


erupts  by  the  reconnection  along  an  arcade  of  magnetic  loops.  In  such  a  case,  the  ends  of 
the  flux  rope  remain  attached  to  the  Sun  during  the  eruption.  A  spheromak  however  is  a 
self  contained  magnetic  field  configuration  supported  by  internal  currents,  with  no  external 
attachment,  and  presumably  if  appropriate,  must  emerge  fully  formed  from  beneath  the  solar 
photosphere. 

In  general,  the  models  in  which  a  flux  rope  forms  by  reconnection  predict  an  acceleration 
phase  at  the  onset  of  flare  reconnection.  By  contrast,  the  2007  May  19  CME  appears 
to  have  moved  out  from  very  low  down  at  almost  constant  speed.  There  is  no  obvious 
feature  in  its  height-time  plot  that  would  correspond  to  the  onset  of  an  acceleration  phase 
coupled  with  reconnection  to  form  a  flux  rope,  suggesting  that  the  observed  flux  rope  must 
have  formed  much  lower  down.  Thus  its  observed  trajectory  suggests  that  a  spheromak 
might  be  a  more  appropriate  description  here,  than  it  would  have  been  for  other  CMEs 
considered  previously  in  iRakowski.  Laming  fc  Lepri  (120071) .  iNakagawa  fc  Matsuokal  (1201011 
find  a  similar  trajectory  for  the  magnetic  cloud  observed  by  ACE  and  Nozomi  on  1999  April 
16-18.  They  also  conclude  that  the  magnetic  cloud  is  better  fitted  by  a  toroidal,  as  opposed 
to  cylindrical,  flux  rope,  though  the  model  fitted  is  still  approximate.  Both  flux  ropes  and 
spheromaks  conserve  hclicity  as  they  expand,  and  to  do_so  must  decrease  their  magnetic 
energy,  a  part  of  which  can  end  up  dissipated  as  heat.  Kumar  fc  Rust  (1199611  consider  an 
expanding  flux  rope,  and  find  temperatures  reaching  1  —  2  x  106K.  We  apply  similar  ideas, 
derived  within  the  context  of  a  spheromak,  to  the  2007  May  19  CME.  We  assume  that  the 
expansion  velocity  of  the  spheromak  entering  the  expression  for  the  heating  is  given  by  the 
observed  motion  of  the  CME. 


5.  Spheromak  non-equilibrium  ionization  model  of  the  19  May  2007  CME 


As  in  IRakowski.  Laming  &  Lepri  (120071)  we  modeled  the  charge  states  within  the  CME 
ejecta  for  a  variety  of  ions  using  an  adaptation  of  the  BLASPHEMER  (BLASt  Propaga¬ 
tion  in  Highly  EMitting  EnviRonment)  code  (ILaming  fc  Grnnll2002U2003l:  lha/ming  fc  Hwang 


2003),  which  follows  the  time  dependent  ionization  balance  and  temperatures  of  a  Lagrangian 
plasma  parcel  as  it  expands  in  the  solar  wind.  The  fundamental  equations  are  outlined 
below.  Further  details  of  the  ionization  and  recombination  calculation  can  be  found  in 
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Rakowski.  Laming  &  Lepri  (120071)  and  references  therein. 

The  density  riiq  of  ions  of  element  i  with  charge  q  is  given  by 
dm, 


nq 


dt 


ne  (Cion^q—iTli  q—i  CionjqTliq^  T  Tle  (Crr^q- |_i  T  Cdr,q+ 1)  ?li  q+1  ^-e  {Crr,q  T  Cdr,q ) 


iq 


(15) 

where  Cion^q,  Crr )(J,  Cdr,q  are  the  rates  for  electron  impact  ionization,  radiative  recombination 
and  dielectronic  recombination  respectively,  out  of  the  charge  state  q.  These  rates  are  the 
same  as  those  used  in  the  recent  ionization  balance  calculations  of  iBrvans  et  al.l  (120061).  with 
more  recent  updates  given  in  Rakowski.  Laming  &  Lepri  (I20071).  The  electron  density  ne  is 
determined  from  the  condition  that  the  plasma  be  electrically  neutral.  The  ion  and  electron 
temperatures,  Tiq  and  Te  are  coupled  by  Coulomb  collisions  by 


dTiq  ( Tiq  -  Te)  q37iiq/  (q  +  1)  ( In  A 

=  —  (J.13n(  1 


dt 


MiqTe'2 


Eiq  7M 


and 


dTe  _  0.13 ne  ^  {Tiq  -  Te)  q2niq/  {q  +  1)  /In  A 


dt 


Tf 


3/2 


E 


Mi 


iq 


iq 


E 


iq 


v  37 


Te  ( dne 
nP 


37 


(16) 


\  dt 


2  dQ 


3nek b  dt  1.5 nPk 


+ 


e^B 

(17) 

Here  Miq  is  the  atomic  mass  of  the  ions  of  element  i  and  charge  q  in  the  plasma,  and  In  A  ~  28 
is  the  Coulomb  logarithm.  The  term  in  dQ/dT  represents  plasma  energy  losses  due  to 
ionization  and  radiation.  Radiation  losses  can  be  taken  from  [Summers  fe  McWhirter  (119791). 
In  a  CME  from  a  filament  eruption  the  densities  and  temperatures  are  such  that  radiative 
losses  can  be  important,  unlike  most  applications  in  the  solar  wind  where  they  are  generally 
negligible.  The  term  —  ( Te/ne )  ( dne/dt)ion  gives  the  reduction  in  electron  temperature  when 
the  electron  density  increases  due  to  ionization.  Recombinations,  which  reduce  the  electron 
density,  do  not  result  in  an  increase  in  the  electron  temperature  in  low  density  plasmas,  since 
the  energy  of  the  recombined  electron  is  radiated  away  (in  either  radiative  or  dielectronic 
recombination),  rather  than  being  shared  with  the  other  plasma  electrons  as  would  be  the 
case  for  three-body  recombination  in  dense  plasmas.  We  include  the  last  electron  heating 
term  to  model  the  Ohmic  dissipation  given  by  equation  10  for  the  heating  in  the  spheromak. 


In  addition  to  ionization  and  recombination  there  is  also  the  geometry  to  consider  in 
the  density  and  temperature  evolution.  For  lack  of  a  better  motivated  profile  we  chose  a 
uniform  initial  density  and  temperature,  exploding  from  an  initial  radius  of  (^o  —  1)  where  zQ 
is  the  starting  height  from  the  center  of  the  sun  in  solar  radii.  Thus  the  underlying  density 
evolution  for  adiabatic  expansion  goes  as  \{z  —  l)/(^o  —  1)]2>  which  goes  over  to  ( z/z0 )2  as 
z  >>  Zq  and  the  CME  expansion  more  closely  follows  that  of  the  ambient  solar  wind. 
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There  are  two  questions  we  are  trying  to  answer  with  our  simulations.  First,  how  well 
can  a  spheromak  solution  explain  the  charge  states  seen  in  the  19  May  2007  event?  To  this 
end,  two  sets  of  models  were  explored.  In  model  1  we  attempted  to  match  both  the  filament 
heating  and  the  ICME  charge  states  using  only  the  heating  dne  to  anomalous  resistivity  in 
the  spheromak.  In  model  2  we  allowed  the  plasma  to  already  have  been  heated  during  the 
eruption  to  temperatures  that  would  explain  the  EUVI  observations,  but  then  followed  the 
subsequent  heating  in  the  spheromak.  Second,  given  that  fast  CMEs,  particularly  at  solar 
maximum,  often  show  Fe  charge  states  as  high  as  Fe1/+  in  situ ,  what  are  the  highest  Fe 
charge  states  obtainable  by  anomalous  resistive  heating  in  a  spheromak  model?  In  other 
words,  how  general  could  such  a  model  be? 

The  baseline  model  1  starts  the  plasma  at  a  temperature  of  104K,  a  density  of  1010cm-3, 
and  a  velocity  of  10  km  s_1  accelerating  to  500  km  s_1.  These  were  chosen  as  appropriate 
temperatures  and  densities  for  cool  filament  material,  and  the  evolution  is  consistent  with 
the  final  density  seen  in  situ.  We  initially  explored  the  effect  on  the  final  charge  states 
of  changing  the  anomalous  resistivity,  and  thus  the  amount  of  magnetic  flux  dissipation, 
through  the  parameter  m.  We  assumed  heating  at  the  average  rate  for  the  spheromak 
volume.  Values  of  m  of  0.01  and  below  leave  the  charge  states  essentially  unaffected  by 
the  spheromak  heating.  (The  charge  states  still  evolve  with  the  adiabatic  expansion  which 
changes  the  temperature  and  density  as  the  CME  erupts).  At  m  values  between  1.5  and 
2  one  reaches  the  limit  of  diminishing  returns  whereby  increasing  m  does  not  result  in 
higher  charge  states,  due  to  the  factor  (R0/R)2"1  in  equation  10.  This  corresponds  to  an 
anomalous  resistivity  between  8  x  10~5s  and  10”4s,  as  much  as  1011  times  the  Spitzer  value, 
with  a  maximum  Fe  charge  state  of  Fe13+.  We  note  that  increasing  the  final  velocity  or 
the  acceleration  rate  does  not  substantially  change  this  limit.  A  more  dramatic  spheromak 
expansion  starting  from  a  smaller  sphere  does  heat  the  material  faster,  but  it  also  reduces  the 
density  too  quickly  such  that  (a)  the  final  densities  are  substantially  lower  than  seen  in  situ 
and  (b)  ionization  declines  such  that  the  highest  charge  states  are  not  produced.  Thus  we 
conclude  that  the  high  charge  states  over  Fe13+  in  fast  CMEs  cannot  be  reproduced  purely 
with  the  spheromak  solution  and  anomalous  resistivity.  This  leaves  the  question  of  how  well 
the  spheromak  can  explain  this  particular  19  May  2007  CME. 

We  choose  m  =  0.75  and  compute  ion  charge  states  averaged  over  the  whole  spheromak 
assuming  the  distribution  of  heating  rates  in  Figure  3.  Since  the  2007  May  22  CME  flux 
rope  does  not  exhibit  Fe  more  ionized  than  Fe13+,  its  highest  charge  states  can  be  explained 
with  the  baseline  model  1  and  an  average  m  =  0.75.  The  overall  weighted  charge  states 
of  this  model  are  shown  in  Figure  5.  There  are  two  obvious  discrepancies  between  in  situ 
observations  and  the  Fe  charge  states  modeled  here.  First,  there  are  too  many  low  charge 
states  from  the  unheated  portion  of  the  spheromak.  This  is  not  a  major  concern  as  it  is 
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possible  that  that  portion  of  the  spheromak  is  simply  not  filled  with  plasma.  Second,  Carbon 
charge  states  seem  to  indicate  that  the  plasma  started  hotter  than  in  these  simulations. 

Additionally  one  can  consider  whether  the  spheromak  heating  could  explain  the  EUV 
observations  of  the  filament  around  launch.  Figure  6  shows  the  early  evolution  of  the  charge 
states  corresponding  to  the  average  heating  rate  near  the  sun  in  the  baseline  model  1  with 
m  =  0.75.  This  started  with  the  plasma  at  l.li2©,  above  the  heights  of  1.07-R©  where 
Liewer  et  al.l  (2009)  first  see  Fe  IX  and  Fe  XII  emission.  Starting  simulations  below  this 


height,  with  a  density  that  evolves  to  become  the  observed  density  in  situ  at  earth  orbit, 
leads  to  strong  radiative  cooling  that  keeps  the  electron  temperature  below  104K  until 
a  radius  of  1.1_R0  is  reached,  whereupon  the  density  is  sufficiently  reduced  to  lengthen  the 
radiative  cooling  time.  Hence  the  heating  observed  by  Liewer  et  al.  (20091).  beginning  before 
the  actual  eruption  began,  must  have  had  a  different  origin. 

Testing  models  that  started  at  lower  heights,  we  found  that  indeed  Fe  IX  and  even  Fe 
XII  could  be  produced  within  a  short  distance  of  launch.  We  did  not  use  these  for  our  in  situ 
modeling  because  for  reasonable  starting  filament  densities  the  final  density  due  to  the  rapid 
expansion  was  too  low.  Overall  we  conclude  that  while  reproducing  the  EUV  observations 
may  be  possible  within  a  spheromak  context,  given  the  sequence  of  heating  and  cooling  of 
the  filament  leading  up  the  event  as  described  by  Bone  et  al.  (2009),  it  is  also  plausible,  and 
probably  more  likely,  that  whatever  mechanism  was  heating  the  filament  pre-eruption  was 
responsible  for  the  filament  heating  during  the  eruption. 

The  above  discussion  naturally  leads  us  to  our  second  model  to  determine  if  spheromak 
heating  starting  from  temperatures  seen  in  EUV  observations  can  explain  the  in  situ  charge 
states.  In  Model  2  we  started  the  plasma  at  106K,  109cm~3,  and  a  velocity  of  100  km  s_1 
similar  to  the  velocity  of  the  filament  from  flLiewer  et  al.l  20091)  allowing  it  to  then  accelerate 
to  the  coasting  velocity  of  500  km  s-1.  Starting  from  this  high  temperature,  the  additional 
heating  of  the  spheromak  has  no  noticeable  effect  on  the  charge  states.  While  the  tempera¬ 
ture  was  chosen  such  that  the  charge  states  at  low  heights  match  observations,  the  final  Fe 
charge  states  fall  short  of  those  detected  in  situ.  We  also  experimented  with  starting  at  106K 
but  the  higher  density  and/or  at  a  slower  velocity.  However,  radiative  recombination  domi¬ 
nated  bringing  the  charge  states  and  temperature  rapidly  back  to  their  104K  levels.  This  has 
implications  for  the  rapidity  and  amount  of  energy  that  had  to  be  dumped  into  the  filaments 
to  produce  their  high  charge  states.  A  thorough  examination  of  this  phenomenology  is  war¬ 
ranted  but  we  leave  this  to  future  work.  Overall,  starting  from  a  higher  initial  temperature 
in  a  uniformly  filled  spheromak  did  not  solve  the  problem  of  matching  the  highest  Fe  charge 
states. 
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6.  Discussion  and  Conclusions 


We  have  explored  modeling  this  filament  eruption  as  due  to  the  heating  provided  in  an 
expanding  spheromak.  Figure  6  shows  the  evolution  of  H,  C,  O,  Si  and  Fe  charge  states 
obtained  with  an  anomalous  resistivity  increased  from  the  classical  (i.e.  Spitzer)  value  by 
a  factor  around  2  x  107  —  6  x  1011,  as  the  plasma  temperature  varies  between  104  K  and 
10'  K.  Anomalous  resistivity  is  thought  to  arise  in  conditions  where  plasma  turbulence  can 
scatter  current  carrying  electrons,  and  can  do  so  more  effectively  than  the  Coulomb  collisions 
between  electrons  and  ions.  It  may  be  estimated  by  replacing  the  isotropization  frequency 
by  Coulomb  collisions  with  the  relevant  wave-particle  interaction  rate. 


Solutions  of  the  one  dimensional  Vlasov  equation  designed  to  model  the  Buneman  insta¬ 
bility  in  a  reconnecting  current  sheet  flWu  &  Huang  2QQ9),  and  therefore  model  electrons  in 
Langmuir  turbulence,  give  quiescent  values  of  the  anomalous  resistivity  of  xlO'  the  Spitzer 
value,  with  transient  values  as  high  as  xlO9  in  a  background  plasma  temperature  of  10'  K. 
In  our  case,  dropping  the  anomalous  resistivity  to  xlO10  the  classical  value,  the  dominant 
Fe  charge  state  becomes  Fe'+.  The  ratio  of  the  electron  plasma  frequency  to  the  collision 
frequency  at  1.4i?0  in  Figure  6  evaluates  to  ~  109,  and  probably  represents  the  maximum 
enhancement  in  resistivity  than  one  may  reasonably  expect.  Lin  et  al.  (120071)  estimate  a 
much  higher  resistivity  from  the  observed  width  of  current  sheets  trailing  CMEs,  as  high 
as  1012  times  the  classical  resistivity.  IBemporad  (12008!)  revisits  this  and  argues  that  what 
is  observed  as  a  macroscopic  current  sheet  is  in  fact  an  assembly  of  microscopic  current 
sheets,  each  with  width  corresponding  to  anomalous  resistivities  similar  to  those  given  by 
(jWu  &  Huang  2009),  although  possibly  attributed  to  different  modes  of  turbulence.  We 
therefore  conclude  that  magnetic  energy  dissipated  in  the  expansion  of  a  spheromak  may  ac¬ 
count  for  the  filament  heating  between  104  K  and  106  K,  but  is  unlikely  to  be  the  explanation 
for  the  high  Fe  charge  states  characteristic  of  temperatures  of  order  10'  K  and  higher.  The 
nature  of  the  heating  specified  in  the  spheromak  solution  is  such  that  the  heating  rate  due 
to  anomalous  resistivity  cannot  be  increased  without  limit,  even  if  an  anomalous  resistivity 
as  high  as  allowed  by  the  solution  could  be  justified  physically. 

Higher  charge  states,  beyond  Fe13+  must  therefore  result  from  other  means  of  heating. 
Noting  that  in  the  2007  May  19  CME,  the  highest  Fe  charge  states  are  observed  exterior  to  the 
flux  ropes,  we  speculate  that  reconnection  associated  with  the  eruption  is  most  likely.  Type 


II  and  Type  III  radio  bursts  were  observed  during  this  eruption  (Kerdraon  et  al 


times  when  the  CME  shock  was  between  1.2  and  1.4  R&  heliocentric  distance  (IVeronig  et  al. 


2010  at 


2008|).  Bearing  in  mind  that  the  plasma  we  observe  and  model  will  be  at  lower  altitude 
that  the  CME  shock,  this  corresponds  very  well  to  the  epochs  of  heating  shown  in  Figure 
6.  Recent  results  suggest  that  reconnection  might  be  an  efficient  means  of  heating  electrons 
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(Drake  et  al. 

2006; 

Oka  et  al. 

2010) 

In  earlier  work  (IRakowski.  Laming  fc  Lepri  1120071).  we  expressed  the  heat  input  to  the 
CME  in  terms  of  its  kinetic  energy  input  by  relating  the  heating  to  the  acceleration,  tac¬ 
itly  assuming  that  the  mechanism  of  acceleration  inevitably  also  supplies  heat  to  the  CME 
plasma.  In  the  case  of  the  spheromak,  the  heating  is  related  to  the  expansion,  but  not  the 
acceleration,  so  the  ratio  of  heat  to  kinetic  energy  does  not  have  a  natural  interpretation. 


We  estimate  it,  though,  to  facilitate  comparison  with  Rakowski.  Earning  &  Lepri 
From  equations  11  and  12, 


(20071). 


/  edt  0.033-Bq.R3  In  (R/R0) 


—  0.2^1  In 

V 


R_ 

Ro 


(18) 


mv2  / 2  2npR3v2/3 

which  close  to  the  Sun  evaluates  to  ~  1.  This  is  similar  to,  but  slightly  lower  than  typical 


values  found  bv  IRakowski.  Laming  fc  Lepri  (120071).  The  CMEs  studied  in  that  work  fre¬ 
quently  had  charge  states  of  Fe  up  to  neon-like,  rather  higher  than  in  the  flux  rope  studied 
here  which  appear  to  be  consistent  with  our  estimate  of  thermal  to  kinetic  energy  in  the 
CME. 


This  work  has  been  supported  by  NASA  Contract  NNG08EK62I  and  by  basic  research 
funds  of  the  Office  of  Naval  Research. 


REFERENCES 

Akmal,  A.,  Raymond,  J.  C.,  Vourlidas,  A.,  Thompson,  B.,  Ciaravclla,  A.,  Ko,  Y.-K.,  Uzzo, 
M.,  &  Wu,  R.  2001,  ApJ,  553,  922 

Antiochos,  S.  K.,  DeVore,  C.  R.,  &  Klimchuk,  J.  A.  1999,  ApJ,  510,  485 
Bemporad,  A.  2008,  ApJ,  689,  572 

Bone,  L.  A.,  van  Driel-Gesztelyi,  L.,  Culhane,  J.  L.,  Aulanier,  G.,  &  Liewer,  P.  2009, 
Sol.  Phys.,  259,  31 

Bryans,  P.,  Badnell,  N.  R.,  Gorczyca,  T.  W.,  Laming,  J.  M.,  Mitthumsiri,  W.,  &  Savin,  D. 
W.  2006,  ApJS,  167,  343 

Ciaravella,  A.,  Raymond,  J.  C.,  Reale,  F.,  Strachan,  L.,  &  Peres,  G.,  2001,  ApJ,  557,  351 
Chandrasekhar,  S.,  &  Kendall,  P.  C.  1957,  ApJ,  126,  457 


Drake,  J.  F.,  Swisdak,  M.,  Che,  H.,  &  Shay,  M.  A.  2006,  Nature,  433,  553 
Fan,  Y.,  &  Gibson,  S.  E.  2004,  ApJ,  609,  1123 
Gloeckler,  G.  et  ah  1998,  Space  Science  Reviews,  86,  497 
Gopalswamy,  N.,  et  ah  2009,  ApJ,  691,  L123 

Henke,  T.,  Woch,  J.,  Schwenn,  R.,  Mall,  U.,  Gloeckler,  G.,  von  Steiger,  R.,  Forsyth,  R.  J., 
&  Balogh,  A..  2001,  J.  Geophys.  Res.,  106,  10597 

Henke,  T.,  Woch,  J.,  Mall,  U.,  Livi,  S.,  Wilken,  B.,  Schwenn,  R.,  Gloeckler,  G.,  von  Steiger, 

R. ,  Forsyth,  R.  J.,  &  Balogh,  A.  1998,  Geophys.  Res.  Lett.,  25,  3465 

Kataoka,  R.,  Ebisuzaki,  T.,  Kusano,  K.,  Shiota,  D.,  Inoue,  S.,  Yamamoto,  T.  T.,  &  Toku- 
maru,  M.  2009,  J.  Geophys.  Res.,  114,  A10102 

Kerdraon,  A.,  Pick,  M.,  Hoang,  S.,  Wang,  Y.-M.,  &  Haggerty,  D.,  2010,  ApJ,  715,  468 

Kilpua,  E.  K.  J.,  et  al.  2009,  Sol.  Phys.,  254,  325, 

Kumar,  A.,  &  Rust,  D.  M.  1996,  J.  Geophys.  Res.,  101,  15667 

Laming,  J.  M.,  &  Grun,  J.  2002,  Phys.  Rev.  Lett.,  89,  125002 

Laming,  J.  M.,  &  Grun,  J.  2003,  Physics  of  Plasmas,  10,  1614 

Laming,  J.  M.,  &  Hwang,  U.  2003,  ApJ,  597,  347 

Lee,  J.-Y.,  Raymond,  J.  C.,  Ko,  Y.-K.,  &  Kim,  K.-S.  2009,  ApJ,  692,  1271 

Lepri,  S.  T.,  &  Zurbuchen,  T.  H.  2004,  J.  Geophys.  Res.,  109,  A1112 

Lepri,  S.  T.,  Zurbuchen,  T.  H.,  Fisk,  L.  A.,  Richardson,  I.  G.,  Cane,  H.  V.,  &  Gloeckler,  G. 
2001,  J.  Geophys.  Res.,  106,  29231 

Li,  Y.,  Lynch,  B.  J.,  Stenborg,  G.,  Luhmann,  J.  G.,  Huttunen,  K.  E.  J.,  Welsch,  B.  T., 
Liewer,  P.  C.,  &  Vourlidas,  A.  2008,  ApJ,  681,  L37 

Liewer,  P.  C.,  de  Jong,  E.  M.,  Hall,  J.  R.,  Howard,  R.  A.,  Thompson,  W.  T.,  Culhane,  J.  L., 
Bone,  L.,  &  van  Dricl-Gesztelyi,  L.  2009,  Sol.  Phys.,  256,  57 

Lin,  J.,  Ko,  Y.-K.,  Sui,  L.,  Raymond,  J.  C.,  Stenborg,  G.  A.,  Jiang,  Y.,  Zhao,  S.,  &  Mancuso, 

S.  2005,  ApJ,  622,  1251 


-16- 

Lin,  J.,  Li,  J.,  Forbes,  T.  G.,  Ko,  Y.-K.,  Raymond,  J.  C.,  &  Vourlidas,  A.  2007,  ApJ,  658, 
L123 

Linker,  J.  A.,  Mikic,  Z.,  Lioncllo,  R.,  Riley,  P.,  Amari,  T.,  &  Odstrcil,  D.  2003,  Phys. 
Plasmas,  10,  1971 

Liu,  Y.,  Luhmann,  J.  G.,  Huttunen,  K.  E.  J.,  Lin,  R.  P.,  Bale,  S.  D.,  Russell,  C.  T.,  & 
Galvin,  A.  B.  2008,  ApJ,  677,  L133 

Lynch,  B.  J.,  Zurbuchen,  T.  H.,  Fisk,  L.  A.,  &  Antiochos,  S.  K.  2003,  J.  Geophys.  Res.,  108, 
A6,  1239 

Lynch,  B.  J.,  Antiochos,  S.  K.,  MacNeice,  P.  J.,  Zurbuchen,  T.  H.,  &  Fisk,  L.  A.  2004,  ApJ, 
617,  589 

Lyutikov,  M.,  &  Gourgouliatos,  K.  N.  2010,  Solar  Physics,  submitted,  arXiv:  1009.3463 
Magara,  T.,  &  Longcope,  D.  W.  2003,  ApJ,  586,  630 

Moore,  R.  L.,  Sterling,  A.  C.,  Hudson,  H.  S.,  &  Lemen,  J.  R.  2001,  ApJ,  552,  833 

Nakagawa,  T.,  &  Matsuoka,  A.  2010,  J.  Geophys.  Res.,  in  press 

Oka,  M.,  Phan,  T.-D.,  Krucker,  S.,  Fujimoto,  M.,  &  Shinohara,  I.  2010,  ApJ,  714,  915 

Rakowski,  C.  E.,  Laming,  J.  M.,  &  Lepri,  S.  T.  2007,  ApJ,  667,  602 

Sui,  L.,  &  Holman,  G.  2003,  ApJ,  596,  L251 

Sui,  L.,  Holman,  G.  D.,  White,  S.  M.,  &  Zhang,  J.  2005,  ApJ,  633,  1175 
Summers,  H.  P.,  &  McWhirter,  R.  W.  P.,  1979  J.  Phys.  B.,  12,  2387 
Veronig,  A.  M.,  Temmer,  M.,  &  Vrsnak,  B.  2008,  ApJ,  681,  L113 
Wu,  G.  P.,  &  Huang,  G.  L.  2009,  A&A,  502,  341 

Zhang,  J.,  Dere,  K.  P.,  Howard,  R.  A.,  Kundu,  M.  R.,  &  White,  S.  M.  2001,  ApJ,  559,  452 
Zhang,  J.,  Dere,  K.  p.,  Howard,  R.  A.,  &  Vourlidas,  A.  2004,  ApJ,  604,  420 
Zhang,  J.,  &  Dere,  K.  P.  2006,  ApJ,  649,  1100 

Zurbuchen,  T.  H.,  Gloeckler,  G.,  Ipavich,  F.,  Raines,  J.,  Smith,  C.  W.,  &  Fisk,  L.  A.  2004, 
Geophys.  Res.  Lett.,  31,  11805 


17  - 

Zurbuchen,  T.  H.,  &  Richardson,  I.  G.  2006,  Space  Science  Reviews,  123,  31 


Fig.  1. —  Fe  charge  state  distributions  detected  by  STEREO  A  and  B  during  the  2007  May 
19  CME  passage.  Highest  Fe  charge  states  are  detected  by  STEREO  A,  though  the  magnetic 
cloud  is  more  clearly  detected  in  STEREO  B.  Magnetic  field  and  other  data  indicate  cloud 
passage  between  3:36  hr  and  16:34  hr  on  day  142  (22  May)  for  STEREO  B  and  21  May  19:12 
until  22  May  00:14  for  STEREO  A  (Liu  et  al.  20081.  The  dashed  vertical  white  lines  show 
the  temporal  extent  of  the  magnetic  cloud  (or  flux  rope)  in  this  event. 
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Fig.  2. —  ACE  charge  state  data  for  C,  O,  Si  and  Fe  in  8  hr  bins.  Earliest  times  have  been 
offset  upwards  for  display  purposes.  The  highest  Fe  charge  states  only  appear  starting  in 
Day  143  after  the  passage  of  the  interior  of  the  magnetic  cloud.  Liu  et  ah  (20081)  report  the 
1CME  passages  through  ACE  as  21  May  22:19  to  22  May  12:43. 


Fig.  3. —  Normalized  dissipation  rate  e  as  function  of  radial  distance  r/R(t )  and  polar  angle 
6  for  two  self-similar  dissipative  structures,  m=l.  Left  Panel :  solutions  ([MS]).  Right  Panel : 
solutions  flZHHD. 
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Fig.  4. —  Distribution  of  volumetric  heating  rates  in  the  spheromak  model  of  the  first  type, 
solution  ([MED,  i11  terms  of  the  average  heating  rate  over  the  whole  volume. 
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Fig.  5. —  Final  charge  state  distributions  of  C,  0,  Si  and  Fe  for  model  1.  The  individual  lines 
show  the  contributions  from  regions  within  the  spheromak  with  different  heating  rates.  The 
histograms  show  the  cumulative  charge  distribution  assuming  a  volume  weighted  summation. 
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